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We show that the cosmic ray (CR) knee can be entirely explained by energy-dependent CR 
leakage from the Milky Way, with an excellent fit to all existing data. We test this hypothesis 
calculating the trajectories of individual CRs in the Galactic magnetic field. We find that the CR 
escape time Tesc{E) exhibits a knee-like structure around E/Z = few x 10^® eV for small coherence 
lengths and strengths of the turbulent magnetic field. The resulting intensities for different groups of 
nuclei are consistent with the ones determined by KASCADE and KASCADE-Grande, using simple 
power-laws as injection spectra. The transition from Galactic to extragalactic CRs is terminated 
at « 2 X 10^® eV, while extragalactic CRs contribute significantly to the subdominant proton flux 
already for > 2 x 10^® eV. The natural source of extragalactic CRs in the intermediate energy region 
up to the ankle are in this model normal and starburst galaxies. The escape model provides a good 
fit to ln(A) data; it predicts that the phase of the CR dipole varies strongly in the energy range 
between 1 x 10^^ and 3 x 10^® eV, while our estimate for the dipole magnitude is consistent with 
observations. 


I. INTRODUCTION 


The all-particle cosmic ray (GR) energy spectrum is 
a nearly featureless power-law between ~ 10^° eV and 
~ 10^° eV, with only a few breaks in its spectral in¬ 
dex. The two most prominent ones are the knee at 
E\^ K, 4PeV, and the ankle at Eg, ~ 4EeV. They must 
contain information about either GR sources or GR prop¬ 
agation. The range of possible theoretical explanations 
for these two breaks has been reduced in the last decade, 
but there is still no firm consensus on their origins. An¬ 
other related open question is where the transition from 
Galactic to extragalactic CRs occurs. Unveiling this 
transition energy holds precious keys to understanding 
particle acceleration in the Universe. Other features have 
been observed in the spectrum between E^ and Eg,, see 
6-g- should contain additional clues. 

In addition to the all-particle spectrum, both measure¬ 
ments of the primary composition and upper limits on 
the CR dipole anisotropy are crucial to constrain models 
of the knee and of the intermediate region up to the an¬ 
kle. The KASCADE-Grande collaboration has recently 
provided measurements of the intensities of individual 
groups of CR nuclei up to U 10^^ eV Qii- Also, the 
IceTop collaboration has presented measurements of the 
average of the logarithmic mass up to 30PeV while 
the Auger collaboration deduced the contribution of in¬ 
dividual CR groups to the total CR flux above 10^^ ® eV 
from studies of the development of air showers [1,0 • In 
the region between the knee and the ankle, upper limits 
on the amplitude of the anisotropy have been reported, at 
roughly the percent leve l, by KASCADE 0, KASCADE- 
Grande 0, and Auger [l^- Below the knee, the dipole 
amplitude has been measured at the ~ 10“® level, by 
notably Sup er-Kamiokande (lOTeV) [Tl|, MILAGRO 
(6TeV) dl], EAS-TOP 100 and 400TeV) [^Ice- 
Cube (20 and 400 TeV) [H and IceTop (2PeV) [l|. 

For the knee, two main explanations currently remain 


possible. First, it may be the signature of the maximum 
energy to which Galactic CR sources can accelerate pro¬ 
tons, see e.g. M- A nearby source could also leave such 
an imprint in the spectrum 17|. Second, the knee could 
be caused by a change in the energy dependence of the 
CR diffusion coefficient and thence confinement time in 
the Galaxy if the CR Larmor radius is the or¬ 

der of the coherence length Ig of the turbulent Galac¬ 
tic magnetic field (GMF) at E\^. In Ref. [2^, we have 
studied this possibility—which we denote as the “escape 
model”—by propagating individual CRs in recent GMF 
models. This enabled us to avoid limitations from the 
diffusion approximation; While reliable analytical ap¬ 
proximations for the diffusion tensor are only available 
in certain limiting cases, the diffusion approximation per 
se is not justified at the highest energies studied. We 
showed that the escape model is viable and can explain 
the individual fluxes of CR groups as measured by KAS¬ 
CADE and KASCADE-Grande. Moreover, our estimate 
for the dipole anisotropy in this model was consistent 
within uncertainties with observations. 


In this work, we extend our previous study in Ref. [2^ 
and formulate a model for the entire energy region be¬ 
tween 300GeV/Z and the ankle. In addition to the 
Jansson-Farrar (JF) GMF model [2lj used in we 
consider the Pshirkov et al. (PTKN) model [2^. 
This enables us to check the dependence of our results 
on the specific GMF model. Moreover, a more detailed 
study of the Galactic CR primary composition between 
Uk and Eg, is presented here and compared to observa¬ 
tions. We show that any remaining heavy nuclei flux in 
the sub-ankle region would be dominated by only one 
or few local sources. We use limits on the iron frac¬ 
tion at > 7 X 10^^ eV determined by the Auger collab¬ 
oration together with ln(A) measurements to constrain 
the transition energy between Galactic and extragalactic 
CRs, deducing 77.niax = £^max,Fe/(26e) ^ 10^^ V as the 
maximal rigidity 77.niax to which Galactic CR sources are 
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able to accelerate CRs. The recovery of the proton and 
helium spectra above EjZ ^ 10^® eV in the KASCADE- 
Grande data is maninly explained by the specific shape of 
the escape rate Tesc(A) discovered in [2^. We show also 
that observational constraints from anisotropy limits are 
compatible with the escape model. A natural extension 
of the ’escape model’ to other normal galaxies suggests 
that the extragalactic flux in the intermediate energy re¬ 
gion up to the ankle is composed of CRs accelerated in 
starburst galaxies. 

This article in organized as follows: In Section m we 
test two GMF models for the regular and turbulent fields, 
as well as different strengths and coherence lengths for 
the turbulence. We deduce a range of models that fit 
constraints from notably the B/C ratio. We then com¬ 
pute in Section uni the resulting fluxes of Galactic CR 
groups and show that they fit very well KASCADE and 
KASCADE-Grande measurements. In Section IlYl we 
discuss the transition from Galactic to extragalactic CRs 
in our model. Finally, we review in Section |V] the con¬ 
straints on and consequences of the escape model, before 
we present our conclusions in Section [VII 


II. GALACTIC MAGNETIC FIELD MODELS, 
PARAMETERS FOR THE TURBULENCE AND 
CR CONFINEMENT IN THE GALAXY 

A. Grammage 


An important constraint on CR propagation models 
comes from ratios of stable primaries and secondaries 
produced by CR interactions on gas in the Galactic disk. 
In particular, the B/C ratio has been recently measured 
by the AMS-02 experiment up to 670GeV/nucleon [^ . 
Above E > 10 GeV/nucleon, these measurements for the 
B/C ratio are consistent with a straight power-law. 

In our previous work m, we used a fit of the gram¬ 
mage performed in Ref. 2^ using the leaky-box formal¬ 
ism and earlier data. In all cases considered in [^ . 
the grammage traversed by CRs at reference energies 
E^jZ = (5 — 15) GeV was found to lie in the range 
(9 — 14) g/cm^. In order to take advantage of the high- 
quality and the large energy range of the B/C data from 
AMS-02, we use now these data to derive the grammage 
traversed by CRs as function of their energies in two sim¬ 
ple models. In the first one, we approximate the fraction 
of the B to C intensities by 


Psp^s 
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where Xi = nipjat are the interaction lengths (in 
gr/cm^), ai the total inelastic cross section, nip the pro¬ 
ton mass and Psp = Csp/atot is the spallation probability 
deduced from the cross sections given in Ref. [2^. In 
the second approximation, we employ a fit function giv¬ 
ing the B/C ratio directly as function of the grammage, 
following the approach in Refs. [13, H^. In Fig. [U the 


grammage derived in Ref. [25j | using earlier data is shown 
as black cross. The grammage deduced from the AMS-02 
data using the first approximation is shown with magenta 
error-bars, while the grammage obtained using the sec¬ 
ond approximation is shown with blue error-bars. Note 
that the error-bars take into account only the statistical 
and systematic errors of the AMS-02 measurement, while 
uncertainties in the cross sections or deficiencies of our 
approximations are not accounted for. The latter can be 
estimated by the differences between the results from the 
two approximations used. 

In order to compare these measured values to those 
predicted in the escape model, we inject N cosmic rays 
at z = 0 in the Galaxy and follow their trajectories Xi(t) 
until they reach the edge of the Galaxy. As radial distri¬ 
bution of the injection points, we use 

n{r) cx (r/i?0)°-’’exp[-3.5(r-i?Q)/i?0] (2) 

with [n] =kpc“^, assuming that the surface density of 
CR sources follows the distribution of supernova rem¬ 
nants in the Galaxy [1^. Here Rq = 8.5 kpc is the 
distance of the Sun to the Galactic center. We em¬ 
ploy n{z) = no exp(—(z/zi/ 2 )^) as model for the gas 
distribution in the Galactic disk, where z is the dis¬ 
tance to the Galactic plane, no = 0.3/cm^ at Rq and 
-21/2 = 0.21 kpc inspired by [s^. We set n = 10“^/cm^ 
as minimum gas density up to the edge of the Milky Way 
at |z| = 10kpc. Then we calculate the average gram¬ 
mage (A) = J dtp{xi{t)) summing up the 

density along the trajectories of individual CRs. Since 
the grammage X[E) oc E~^ scales as the confinement 
time t{E) oc E~^ , this quantity serves also as an indica¬ 
tor for changes in the CR intensity induced by a variation 
of the CR leakage rate. 


B. Jansson-Farrar model for the GMF 

Let us recall first how the properties of a turbulent 
magnetic field determine the propagation of charged par¬ 
ticles in the diffusion picture, before we discuss the spe¬ 
cific case of the JF model. A turbulent magnetic field is 
characterized by its power-spectrum 'P(fe). The maximal 
length Lniax of the fluctuations and the correlation length 
Ic are connected by = {a — l)Lmax/(2Q;) for V{k) oc 
A:““. Assuming that the turbulent field is isotropic, the 
slope of the power-spectrum V[k) oc k~°‘ determines the 
energy dependence of the diffusion coefficient D{E) in the 
limit E <C Ect as D{E) oc on distances I :s> Lmax. 

Here, the critical energy Ecr is defined by Ri,{E„) = Ic 
and thus the condition E <C Ecr ensures large-angle scat¬ 
tering, while the requirement I ^ Lmax guarantees that 
features of anisotropic diffusion are washed out. Finally, 
we recall that the confinement time r scales as the inverse 
of the diffusion coefficient. 

In our previous work (^ . we used the JF model for 
the regular and turbulent components of the Galactic 
magnetic held [2l|, choosing as the maximal length of 
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the fluctuations Lmax = 10 pc. Note that for Kolmogorov 
turbulence the maximal length of the fluctuations Lmax 
and the correlation length Ic are connected by Lmax = ^Ic 
and that the diffusion coefficient scales as D{E) oc 
for E <C Eci- We considered two values of its root mean 
square (rms) strength, the original one suggested in [2l| 
(/3 = 1) and a second one rescaling it to one tenth of its 
original value {f] = 1/10). 

In Fig. [U we compare the grammage calculated from 
simulated CR trajectories for these two cases with the 
grammage deduced from B/C measurements. Because 
of the large energy reach of the AMS-02 data, the ex¬ 
trapolation required from the lowest energy of our nu¬ 
merical calculations, E = 10^^ eV, to the measurements 
has decreased to two orders of magnitude. Using the JF 
model with /3 = 1 as proposed in [2l| would require a 
constant power spectrum of magnetic field fluctuations, 
V{k) oc k~°‘ with a = 0, in the intermediate energy 
range. Such a power-spectrum is difficult to reconcile 
with the theoretical understanding of turbulence. More¬ 
over, the CR spectrum is very close to a power-law above 
~ 300 GV. This implies that if D{E) would become sig¬ 
nificantly flatter beyond TeV energies (e.g. changing from 
D{E) oc to oc i?°), then the injection spectrum of 
sources has to have the exact opposite change of slope 
(e.g. respectively from oc E~‘^ '^ to oc E~‘^-'^). Alterna¬ 
tively, a change in the source density should compensate 
the change in D{E) such that the observed CR intensity 
remains a nearly featureless power-law [sij. Although 
such a conspiracy cannot be excluded, it appears to us 
as a not very attractive option. 

Choosing a Kolmogorov^ power-spectrum 'P{k) oc 
k~^^^ as the theoretical model with the smallest slope 
a, we have to reduce Rrms therefore by a scaling factor 
/3 < 1 relative to the Rrms suggested in [2l[. The exact 
value of (3 depends on the chosen coherence length 1^: 
A smaller coherence length leads to faster diffusion and 
thus to a smaller value of the grammage. For instance, 
a coherence length close to the upper limits derived in 
Ref. [13, Zc = 5 pc allows with j5 = 1/8 a somewhat 
weaker reduction in the level of the turbulence, cf. the 
blue line in Fig. [T] Increasing the coherence length even 
further to Ic = 30 pc, the scaling factor can be reduced 
to /3 = 1/5, cf. Fig. [H 

Next we examine how the shape of the grammage X as 
function of energy E/Z depends on the two parameters 
Ic and /3. In [13, we discovered a specific shape of X{E) 
that leads not only to a knee-like feature but reproduced 
also the recovery of the proton and helium spectra above 
E/Z 10^® eV, visible in the KASCADE-Grande data. 
From the examples in Figs. [T] and [H it is clear that a too 
strong turbulent field, (3^1, results in knee-like feature 


^ Note that a Kraichnan power-spectrum (a = 3/2) would require 
a stronger rescaling of the turbulent field, leading to a potential 
conflict with synchrotron data. 



E/Z in eV 


FIG. 1: Grammage for different coherence lengths Ic and tur¬ 
bulent fields: red squares Lmax = 10 pc and /3 = 1, black dots 
Umax = 10 pc and P = 0.1, and blue triangles Lmax = 25 pc 
and P = 0.125; all cases for the JF GMF model [13 • Ad¬ 
ditionally we show the grammage deduced from B/C data. 


at too high energy. Compensating a relatively strong 
turbulent field by decreasing the coherence length tapers 
off both the knee-like feature and the recovery, as shown 
by the case Ic = 30 pc in Fig. [2l As a consequence, the 
allowed range of turbulent field strengths and coherence 
lengths is correlated and very restricted, Zc ~ (1 — 10) pc 
and P ~ 1/10 - 1/8. 

This behavior is best illustrated comparing the mod¬ 
ulation induced by the energy dependence of X{E) on 
the intensity of protons to KASCADE and KASCADE- 
Grande data, and asking that a certain parameter choice 
reproduces the position and the shape of the proton 
knee. In Fig. [31 we show I{E) = Io{E)X{E), where in 
Io{E) = Io{E/Eo)°' the normalisation Iq and the slope 
a are chosen such to obtain a good agreement with ob¬ 
servations below the knee for the case of full (/3 = 1, 
Umax = 10 pc) and reduced (/3 = 0.1) turbulent fields 
with Umax = 10 pc: This comparison demonstrates that 
only the case with reduced turbulence can reproduce the 
observed shape of the proton flux. 


C. Pshirkov et al. model for the GMF 

In order to test the dependence of our results on the 
GMF model, we compute additionally the grammage 
X{E) in the Pshirkov et al. model [^J, [l^l- Its regu¬ 
lar field consists of toroidal components in the Galactic 
halo, and of a disk field which follows spiral arms in the 
Galactic plane. We choose the bisymmetric (BSS) bench¬ 
mark model of [13 , where the disk field presents reversals 
between consecutive arms. The authors of [H, [13 did 
not present a three-dimensional model for the turbulent 
field. However, Pshirkov et al. [13 derived upper bounds 
on the deflections of ultra-high energy CR (UHECR) in¬ 
duced by the turbulence magnetic field. Requiring that 
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FIG. 2: Grammage for different coherence lengths Ic and tur¬ 
bulent fields: red squares Lmax = 25 pc and P = 1/8, blue dots 
imax = 150 pc and /? = 1/5, and green triangles Lmax = 25 pc 
and P = 1/5; all cases for the JF model [SJ]. Additionally we 
show the grammage for the PTKN model with Lmax = 25 pc 
and /3 = 0.1 by black stars. 



E[eV] 


FIG. 3: Intensity using I{E) oc X{E) (red solid line) for the 
reduced turbulent field compared to the case of full turbulent 
field (red dashed line). 



2 1.6 1.2 0.8 0.4 0 

FIG. 4: Deflections of 40 EeV protons in a turbulent 
GMF realization with I/max = 100 pc and firms = 

5 fj.Gx exp( —I 2 I/I .8 kpc). 

on the sky, and present them in Fig. 21 after smoothing 
over 5° circles. One can see that 40 EeV proton deflec¬ 
tions in such a turbulence are compatible with the fit for 
the upper limit on deflections from (2^ , both at high and 
low Galactic latitudes. For such a profile, the constraints 
at high latitude are more stringent than those close to 
the Galactic plane. Therefore, the results of Ref. 
imply that, for a scale height of 1.8 kpc of the turbulent 
field, Bq should not be significantly larger than Ri 5^G. 
We take this value in the following. 

We can now compute the grammage in the PTKN 
model for the regular field, which we supplement by our 
toy model for the turbulent field. We find that we have to 
reduce the normalisation Bq of the turbulent field by a si- 
miliar factor P as in the JF model: The case /? = 1/10 and 
Zc = 5 pc is shown in Fig.[2]with black squares. Compared 
to our favorite cases in the JF model {{P = 1/8, Zc = 5pc} 
and {P = 1/10, Zc = 2pc}), the qualitative behavior of 
X{E) is very similar. It is therefore possible to limit the 
numbers of models, and we then decide to use the JF 
model for the rest of this study. 


III. FLUXES OF GALACTIC CR GROUPS 


these bounds are satisfied allows us to construct a toy 
model of the turbulent GMF. 

We choose Zc = 5pc and as profile function firms( 2 ) = 
fioexp(—|z|/1.8kpc). Since the UHECR deflections at 
Earth as predicted by do not show any significant 
dependence on the Galactic longitude, we neglect for sim¬ 
plicity a possible weak dependence of firms on the Galac- 
tocentric radius. 

We backtrack individual 40 EeV protons from the 
Earth in a realization of isotropic Kolmogorov turbulence 
with such characteristics^. We compute their deflections 


A. Diffuse fluxes from all sources 

In order to calculate the GR flux at Earth, we adopt 
the following procedure. First, we compute the three- 
dimensional GR density around a source at different 
times, and later use it as a template to infer the GR 
distributions around other sources. This allows us to 
reduce the required computing time and makes the prob¬ 
lem tractable. For this template, we take a source lo¬ 
cated in the vicinity of the Earth and propagate individ¬ 
ual CR protons from it in the GMF models which are 
presented in the previous Section. We do this for CRs 


^ In order to be compatible with the asumptions of l23ll . we have 


used Ic = 50 pc for this comparison. 
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FIG. 5: Fluxes of CR protons (upper left), helium (upper right), CNO (lower left) and heavy elements including Mg, Si and 


Fe (lower right) are shown in red. Errorbars show the variations in time of the fluxes. Experimental data from PAMELA for 
protons and helium, CREAM [sj, KASCADE [| and KASCADE-Grande 0 ]. 


with energies between 100 TeV and 100 PeV, and take 
four energies per decade. We divide the space around 
the source in cylindrical sectors with radii ranging from 
0.1 kpc to 4kpc. We checked that the remaining contri¬ 
bution from the radial range 4 — 10 kpc does not change 
our results. Their widths in vertical direction depend on 
2 and we choose their boundary lines as \z\ = 100pc, 
200 pc, 300 pc, 400 pc, 500 pc, 700 pc, Ikpc, 2 kpc, 3 kpc, 
5 kpc, and 10 kpc. We divide time in bins of 5kyr with 
a total of 6000 bins, up to t = 30 Myr. We checked also 
that the chosen bin sizes are small enough, in order to 
have no impact on the results. For every propagated 
CR, we save the fraction of its path in a given bin and 
average over all simulated particles. From this, we de¬ 
duce the three-dimensional time-dependent CR density 
in the Galaxy for this source, and then for a distribution 
of sources. 

We create the ensemble of CR sources as follows. We 
generate their positions within the Galactic disk (|z| < 
lOOpc), assuming that the density of GR sources follows 
the distribution given by Eq. This distribution is as¬ 


sumed to depend only on the distance to the Galactic cen¬ 
ter, and not on the direction around it in the plane. The 
remaining parameters are the frequency of GR sources 
and the energy released in CRs by each of them. Only 
the product of these parameters is constrained by ob¬ 
servations. We fix the energy released in CRs by each 
source to Etot = 10®° erg, leaving the source frequency 
as the only free parameter. We sum up the contributions 
from these generated sources to the CR flux at Earth, in 
any time bin, and for a total duration of 300 Myr. For 
each energy bin, we save both the average flux at Earth 
and its one sigma deviation in time, i.e. the values of 
the flux where 16% of the cases are above and below the 
average. Note, that the upper and lower limits are not 
symmetric at high energies, since there are more cases 
with lower flux than average. 

Results for nuclei with charge Z are deduced from the 
above calculations for protons by shifting the energy by 
a factor Z. We then interpolate the resulting CR nuclei 
fluxes to the same energies as for protons. At energies 
below Z X 100 TeV, we assume that diffusion in the Kol- 
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mogorov turbulence shifts the injection power-law oc E~°‘ 
by 1/3 to the spectrum oc observed. This is in¬ 

deed what we observe in our simulations in the energy 
range cii Z x (100 — 300) TeV. Therefore we assume that 
the CR spectrum of protons released by sources follows 
a power-law spectra oc E~^ '^; the maximal energy Ep 
of protons will be fixed later by considering constraints 
form the resulting dipole anisotropy and the observed 
nuclear composition of CRs. For all other nuclei, we use 
power-law spectra with either oc E~^-^'^ or oc and 

maximal energy ZEp. These power-law indices are cho¬ 
sen so as to ht the direct observations from CREAM at 
low energy. We fix the density of sources by normalizing 
the flux found in our simulations to the observed one at 
100 TeV. On average, we require 440 sources per 100 kyr 
for a total energy per source of Etot = 10^° erg, so as to 
fit the observed CR spectra. Within one time bin of 5 
kyr we generate sources according to a Poisson distribu¬ 
tion with an average of 22 sources, corresponding to the 
required source density on larger time scales. 

In Fig. m we plot the CR nuclei fluxes, multiplied by 
as a function of energy. In the upper left and 
upper right panels of Fig. [5l we show the proton and 
helium fluxes, both for turbulent fields with Lmax = 
10 pc and with L„iax = 25 pc. We plot orange lines 
oc £"“ 2 .4-1/3 (upp 0 j. left panel) and oc £- 2 . 22 - 1/3 (upper 
right panel), which represent the slopes expected theoret¬ 
ically at Earth, for our injection spectra with a = 2.4 and 
2.22 as power indices and Kolmogorov turbulence. Note 
that the slopes of the injection spectra required for nu¬ 
clei, (X cx 2.2, coincide with the naive expectations from 
diffusive shock acceleration. Only the proton injection 
spectra requires a somewhat softer slope, a = 2.4, than 
expected. 

In the two upper panels, we show the experimental 
data from PAMELA [s^ grange points), CREAM 
(magenta^ KASCADE Q (green) and KASCADE- 
Grande Q (blue). The proton flux reported by 
KASCADE-Grande is 40% larger than the flux from 
KASCADE in the (10 — 30) PeV region, where the er- 
rorbars of both experiments are relatively small. In con¬ 
trast, the helium flux from KASCADE-Grande is below 
the one measured by KASCADE . This behavior may 
be explained by the insufficient discrimination power be¬ 
tween protons and helium in the KASCADE-Grande ex¬ 
periment [ 3 ^ . For this study, we choose therefore to 
reduce the proton flux of KASCADE-Grande by 40%, 
and add this difference to the helium flux, in same en¬ 
ergy bins. By doing so, the CR fluxes of KASCADE- 
Grande and KASCADE experiments become consistent 
with each other. 

In the lower left panel of Fig. [SJ we plot the CNO flux, 
which predominantly consists of carbon and oxygen. We 
calculate the carbon and oxygen fluxes by normalizing 
them to the CREAM fluxes interpolated to higher en¬ 
ergies with power-laws, and then sum them up. The 
CREAM flux in this hgure is the sum of its carbon and 
oxygen fluxes, where we use carbon energy bins for the 


binning, and interpolate the oxygen flux to these bins be¬ 
fore summing up. KASCADE and KASCADE-Grande 
measurements of the GNO flux are directly compared to 
our fluxes. 

In the lower right panel of Fig. [SJ we show the flux of 
heavy nuclei, which is dominated by Mg, Si and Fe nuclei. 
Since the KASCADE and KASCADE-Grande collabora¬ 
tions divide this flux in two parts ({Mg -|- Si} and Fe), 
we also plot these two contributions separately in Fig. [51 
Mg and Si in the left panel, and Fe in the right panel. 
It is possible to link the last points from CREAM to 
the first points from KASCADE-Grande with a smooth 
power-law, but not to those from KASCADE. This is 
likely to be due to the difficulty for KASCADE to distin¬ 
guish between Si and Fe nuclei . Therefore we choose 
to sum up the Mg, Si and Fe fluxes into a single ’heavy 
component’ in the lower right panel of Fig. [51 In this fig¬ 
ure, the combined ’heavy nuclei flux’ as measured by the 
KASCADE experiment is smooth and agrees well with 
a simple power-law extrapolation of the CREAM flux 
to higher energies. It also agrees with the KASCADE- 
Grande flux. As for the other components, the model 
presented in this work hts well the heavy nuclei flux too. 

As can be seen in Fig. [5] (upper left), the CR pro¬ 
ton flux follows a power-law from 300 GeV up to about 
1 PeV. It then changes to a steeper slope at the knee, 
and recovers at ~ 10 PeV to a flatter power-law with 
index a ~ 2.5. Similar ’knee-like’ cutoffs, shifted by fac¬ 
tors Z in energy, are visible in the fluxes of all groups of 
GR nuclei—see the other panels of Fig. jS] These plots 
demonstrate that the “escape model” hts very well all 
these observations. As discussed previously in the 
knee is due, in this model, to a change in behaviour with 
energy of the CR diffusion coefficient. The energy of the 
knee corresponds to the energy at which the Larmor ra¬ 
dius of CR protons is of the order of the coherence length 
of the turbulent magnetic held (4 = Lmax/5 for a Kol¬ 
mogorov spectrum). For the held strengths we consider 
in this paper, £rms — 0.3/rG (or /3 = 1/8) averaged over 
a circle of radius 1 kpc distance from the solar position in 
the Galactic plane, we hnd in our calculations a change 
in the slope of the CR hux at about 1 PeV, as observed 
in the proton data. 


B. Flux from nearby sources 

In the escape model, the hatter part of the CR pro¬ 
ton hux above ~ 10 PeV is dominated by recent nearby 
sources. This is due to the fact that the conhnement 
time of CR protons in the Galaxy quickly drops with en¬ 
ergy beyond the energy of the knee. The high energy 
part, £ > 3 X 10^^ eV, of the Galactic hux is dominated 
by heavy elements (Mg-|-Si-|-Fe). Recent nearby sources 
would dominate the hux of heavy elements at these en¬ 
ergies. The composition study published by Auger in Q 
constrains however the fraction of iron. Using conserva¬ 
tively the results obtained using the EPOS-LHC simu- 
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FIG. 6: Fluxes of Mg and Si (left) and Fe (right) are shown in red for Zmax = 25 pc and in blue for Zmax = 10 pc. Errorbars 
show the variations in time. Experimental data from CREAM, KASCADE and KASCADE-Grande. 


lation, the iron fraction above 6 x 10^^ eV is limited as 
< 20%. We can add this constraint, excluding all time 
bins where the iron fraction exceeds this Auger limit. In 
the left panel of Fig. [71 we show experimental data for 
iron from CREAM, KASCADE and KASCADE-Grande 
together with the predicted iron flux without (blue error- 
bars) and with accounting of the Auger iron constraint 
(red errorbars). The maximally allowed iron flux using 
the composition constraint from Q is show as a magenta 
line. Since the signature of nearby, recent CR sources 
is a large iron fraction, the Auger constraint effectively 
eliminates these cases, resulting in a reduced flux at high 
energies, E > 10^^ eV. In the right panel of Fig. [71 we 
show the same plot for protons. 

Next we consider the effect of the Auger limit on the 
flux of the source which gives the maximal contribution 
at highest energy 100 PeV. The flux from this dominant 
source is shown in both panels of Fig. [7] without (blue 
circles) and with (red circles) accounting for Auger con¬ 
straint on the iron fraction. In the former case, the dom¬ 
inant source contributes almost 100% of the proton flux 
at 5 X 10^® eV, while taking the Auger iron constraint into 
account reduces the dominance of the strongest source. 
Clearly, the relatively small fraction of iron observed by 
Auger disfavors the presence of a dominating source even 
at the end of Galactic CR spectrum. 

Finally, we note that in the cases when the Auger iron 
constraint is violated, the total proton flux exceeds then 
the measured one and the knee-like structure is less pro¬ 
nounced than in the observed data. Above Z x 10^® eV, 
the predicted CR flux is dominated by a single nearby 
and recent source. Such a situation contradicts not only 
the Auger iron limit, but would violate also the limits on 
the dipole anistropy of the CR flux: If one assumes that 
Galactic sources are able to accelerate only to energy just 
below the Auger constraint on iron, Emax < 7 x 10^^ eV, 
the constribution of recent nearby sources is still limited 
by the Auger anisotropy limits, cf. the discussion in the 


next section. Thus the Auger limits on anisotropy and 
the iron fraction exclude the possibility that the highest 
energy part of the Galactic CR spectrum is dominated 
by a recent nearby source. 


IV. TRANSITION FROM GALACTIC TO 
EXTRAGALACTIC CRS 

Determining at which energy E, the CR flux starts to 
be dominated by extragalactic sources is one of the most 
important unsolved problems in CR physics. While in 
the ’dip model’ of Ref. the transition energy is as 
low as A* Ri a few x 10^^ eV, the ankle has been in other 
models identified with the transition between Galactic 
and extragalactic GRs [13, A* w Aa ri a few x 10^® eV. 
In our model, the energy of the transition A, depends 
both on the maximum rigidity 7?.niax to which Galactic 
sources are able to accelerate CRs and on the distance 
to the nearest active source. As we have seen in the pre¬ 
vious section, the Auger constraint on the iron fraction 
effectively eliminates the possibility that a single source 
dominates the high-energy part of the Galactic CR spec¬ 
trum, reducing thereby the fluctuations. 

We impose in the following the iron constraint through¬ 
out and assume that the end of the Galactic CR spectrum 
is determined by the maximum rigidity T^max to which 
Galactic sources are able to accelerate CRs. Although 
T^max is a free parameter in our model, the energy of the 
transition A* can be determined by using additional ob¬ 
servational constraints. One possibility is to constrain 
the maximum contribution of Galactic sources to the 
total GR flux by using the observational limits on the 
anisotropy of the CR flux. Another way to determine 
A*, is to study the elemental composition of primary CRs 
and to use the fact that the composition of Galactic and 
extragalactic CRs should in principle differ from one an¬ 
other. 
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FIG. 7: Left panel: Experimental data for iron together with the predicted iron flux without (blue errorbars) and with the 
Auger constraint (red errorbars). Additionally, the iron flux from the dominant source is shown without (blue circles) and with 
Auger constraint (red circles). Right panel: The same plot for protons. All cases for /max = 25pc. 



FIG. 8: CR flux and the extragalactic component we pre¬ 
dict for TZmax = 1 X 10^^ eV, shown with the data from 
CREAM 0, KASCADE [sl], KASCADE-Grande 0], TI¬ 
BET [^, and Auger experiments. 


We start with the latter method. As a first step, 
we derive the all-particle CR flux summing up all CR 
groups and compare it to the experimental data of KAS¬ 
CADE KASCADE-Grande [I|, TIBET [H, and 
Auger [^. Then we deduce the extragalactic flux for 
a given TZmax by subtracting the predicted total Galac¬ 
tic flux from the measured total GR flux. The resulting 
extragalactic flux is shown in Fig. [5] with a red solid line 
for TZmax = 1-0 X 10^^ V. Next, we have to fix the nu¬ 
clear composition of the extragalactic CR flux. As a first 
approximation, we can assume that its composition is 
constant in a sufficient small energy interval around E*. 
In contrast, the Galactic GR composition is strongly en¬ 
ergy dependent between the knee and the cutoff of the 
Galactic flux. Thus, we expect that an observable like the 
average of the logarithmic mass number, ln(A), will be 



1e+15 le-fie 1e+17 1e+18 164-19 le-i-20 

E/eV 

FIG. 9: Average of the logarithmic mass In A predicted by 
our model for three different assumptions on the composition 
of extragalactic CRs, versus the experimental data. 


quickly changing for E < while being approximately 
constant for energies slightly above i?*. 

In Fig. [5J we plot measurements of ln(A) from sev¬ 
eral experiments, together with the values of ln(A) calcu¬ 
lated within the ’escape model’ studied here. The points 
for KASGADE have been computed by converting the 
flux measurements given in into ln(A.) values^. The 
most striking feature, namely the peak in ln(A.) around 


® Recall that the KASCADE data for the heavy components 
showed a discrepancy to the extrapolation of the CREAM and 
KASCADE-Grande data and we had to sum them into a sin¬ 
gle heavy component. For the calculation of ln(A), we used in¬ 
stead the original fluxes for the separate CR groups what explains 
the small off-set between our prediction and ln(A) deduced from 
KASCADE data at low energies. 
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5 X 10^® eV, is clearly visible in all data sets, although its 
exact position and strength depend on the experiment. 
Our model reproduces the trend in the data very well. At 
higher energies, the composition becomes lighter because 
of the ’flattening’ of the escape time at such energies, see 
Section im For the value of TZmax we consider here, ex- 
tragalactic CRs start to contribute to the observed flux 
at ss 10^”^ eV. Consequently, above this energy, the ex¬ 
act value and shape of ln(A) depends on the assumed 
composition of the extragalactic flux. In blue, we show 
ln(A) for an extragalactic flux made of protons only, in 
magenta for a mix of 50% p and 50% Fe, and in red for a 
mix of 60% p, 25% He, and 15% N. Independently of the 
composition chosen for the extragalactic component, we 
can identify the energy where ln(A) stops to decrease 
with the maximum energy ifmax.Fe to which Galactic 
sources can accelerate iron. It corresponds to the rigidity 
= ^'max,Fe/26e. This transition is clearly visible in 
the PAO data, around the ankle, and allows us to deter¬ 
mine the maximum rigidity as TZmax — lx 10^^ V. 

Accelerating Galactic CRs to T^max — lx 10^^ V 
is challenging, even for supernovae exploding in dense 
winds. However, several models in which one may reach 
such a high energy have been proposed. One possibility 
is the two-step acceleration of CRs in OB regions, see 
Ref. [dlj for a recent review. Since most of core-collapse 
SNe are located in superbubbles, CRs accelerated by in¬ 
dividual SN remnants may be additionally accelerated 
in superbubbles to energies T^max — lx 10^^ V [dlj. As 
another possibility, Ref. suggests that CRs can be ac¬ 
celerated to ultra-high energies at the termination shock 
of young pulsar winds. Note that in the early stages 
when the acceleration is most effective pulsars stay in the 
same OB regions, and the argument discussed applies in 
this case as well. The TeV gamma-ray emission from ex¬ 
tended Galactic sources was studied in Ref. . There it 
was found that the number of extended sources detected 
in Fermi data is consistent with the expected number of 
TeV GR sources. The majority of these TeV gamma-ray 
sources was associated with pulsars. If these gamma-rays 
have a hadronic origin, pulsars may be candidates for the 
Galactic CR sources. 

For the case of a mixture of 60% p, 25% He, and 15% 
N (red curve in Fig.[^, we obtain a good agreement with 
the ln(A) data from PAO up to 2 x 10^® eV. While this 
choice of composition is not unique, it is consistent with 
the results from the recent composition study published 
in [ 3 ]. In particular. Ref. [ 3 ] found the fraction of iron to 
be below 20% above 6 x 10^^ eV. This agrees well with 
the results of our model, where the Galactic flux at 6 x 
10^"^ eV consists purely of iron but contributes to only 
15% of the total CR flux. 

In addition to fitting the above observables, we still 
have to verify that the model presented here is also 
consistent with the existing upper limits on the CR 
anisotropy. In the diffusion approximation, the CR 
dipole anisotropy d is given by d = 3D Vln(n)/c. Fol¬ 
lowing the same procedure as in [^ . we compute the 



FIG. 10: Dipole amplitude d(E) as a function of energy E 
in the GMF model of Ref. [^ . using a reduced turbulent 
magnetic field with /3 = 1/8 and Lmax = 25 pc. 


average anisotropy and derive the energy dependence of 
D{E) from the escape rate as calculated previously, set¬ 
ting D{E/Z) oc l/Tesc{E/Z). We fix the proportion¬ 
ality constant by requiring that the dipole amplitude 
d = fkdk equals the dipole component d observe d by 
the EAS-TOP collaboration at if = 1.1 x 10^^ eV [l3,l3|. 
Here, k labels the groups of nuclei we consider in the 
Galactic flux plus an extragalactic component. The lat¬ 
ter has a dipole amplitude which is independent of its 
composition and which we set equal to 0.6%, as expected 
for the extragalactic Compton-Getting effect [1^. The 
factor /fe corresponds to the fraction the component k 
contributes to the total CR flux, and dk oc l/Tesc{E/Z) 
to their individual dipole. The relatively low value of the 
CR dipole measurements at TeV-PeV energies is known 
as the ’CR anisotropy problem’. Some authors have sug¬ 
gested that conditions of the local interstellar turbulence 
may be the cause [i^ 113 ] . 

In Fig. [TUI we show the resulting dipole amplitude 
d as a function of energy E. As expected, the ampli¬ 
tude rises below the knee as while it increases 

approximately as until 1 x 10^^ eV. At higher en¬ 
ergies, the dipole amplitude decreases, which is due to 
the facts that the Galactic composition becomes heav¬ 
ier and that the extragalactic contribution grows. We 
also plot the values of d observed by IceGube 0 , as 
well as the 99% G.L . up per limits on d± from the Pierre 
Auger Observatory |10| . Gomparing our estimate for the 
dipole amplitude with the upper limits in the energy 
range 10^^ — 10^® eV, we should take into account that 
the approximation d oc l/resdE/Z) starts to break down 
above E/Z > 10^^ eV, which leads to a sizeable error. 
We conclude therefore that our prediction is marginally 
consistent with these limits. The Pierre Auger Obser¬ 
vatory should however be able to reach a detection of 
the dipole anisotropy. Let us also note that the escape 
model predicts that the phase of the dipole amplitude 
varies strongly in the energy range between 1 x 10^^ and 
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FIG. 11: The extragalactic proton flux deduced from the 
Auger and KASCADE-Grande data is shown in magenta to¬ 
gether with the proton flux observed by KASCADE (green 
errorbars) and KASCADE-Grande (blue errorbars). The to¬ 
tal CR flux from Auger (blue errorbars) is shown together 
with the deduced proton flux (magenta errorbars). 


at higher energies. We show therefore in Fig. [TT] ad¬ 
ditionally the total CR flux (blue) measured by Auger. 
Applying the proton fraction from Ref. Q obtained us¬ 
ing the EPOS-LHC simulation we can then derive the 
resulting proton flux which is shown in magenta. To 
guide the eye, we plotted also a power-law with 

an exponential cutoff aX E = 2.4 x lO'^® eV as an or¬ 
ange line. Such a E~‘^ ‘^ power-law interpolates nicely 
between our prediction for the extragalactic proton flux 
using the KASCADE-Grande data and using the Auger 
data. This conclusion would not change using the spec¬ 
trum and the composition measured by the Telescope Ar¬ 
ray (TA), since the two experiments agree on the points 
that are the most important for our analysis: A very 
small fraction of iron and a large fraction of protons be¬ 
low eV. We conclude therefore that the extragalac¬ 
tic proton flux determined in the escape model, although 
with large errors, is consistent with the slope expected 
from shock acceleration and fits to the proton flux deter¬ 
mined by Auger below the ankle. 


V. DISCUSSION 


3 X lO'^® eV: This corresponds to the range where the tran¬ 
sition from Galactic to extragalactic CRs lies. Such a 
picture is supported by current observations of the phase 
of the dipole, see Refs. EE [ 3 . 

In summary, there are two reasons for having an 
early transition, from predominantly Galactic to predom¬ 
inantly extragalactic CRs, at if a few x lO'^^eV. Eirst, 
the limits on the observed dipole anisotropy requires ei¬ 
ther a very heavy Galactic composition or a predomi¬ 
nantly extragalactic contribution at if > lO'^® eV 3.1^. 
The former possibility is however strongly disfavored by 
the recent composition measurements from the Auger 
collaboration @,0. Second, identifying the energy where 
ln(A) stops decreasing with the maximum energy to 
which Galactic sources can accelerate iron, ifmax.Fe ~ 
3 X I0'^®eV, suggests that the maximal rigidity reached 
in Galactic sources satisfies T^max = ifmax,Fe/(26e) ^ 

Einally, we comment on the contribution of extragalac¬ 
tic protons to the observed proton flux by KASGADE 
and KASGADE-Grande. In Fig.[TTJ we show these exper¬ 
imental data together with the predicted Galactic proton 
flux (red errorbars), taking into account the Auger iron 
constraint. At if > 3 x eV, the predicted Galactic 
proton flux lies below the measured one: Within the es¬ 
cape model, this difference should be accounted for by 
extragalactic protons. Subtracting the measured proton 
flux from the flux calculated in the escape model, we ob¬ 
tain a prediction for the extragalactic proton flux shown 
in magenta in Fig. [TTJ Note that the absolute value of 
this extragalactic proton flux is too small to impact the 
ln(A) plot. Fig. ini We can check if the interpretation 
of this Galactic proton deficit in our model as an extra¬ 
galactic flux makes sense comparing it to expectations 


Before we conclude, we review the main properties of 
the proposed escape model for Galactic CRs and the re¬ 
sulting consequences for the transition between Galactic 
and extragalactic CRs. 


A. Constraints on the GMF 

The “escape model” which aims at explaining the CR 
data from EjZ ^ 300 GeV to 100 PeV by the energy- 
dependent GR leakage from the Milky Way is strongly 
constrained by experimental data: 

1. The position of the knee, ifk « 4 PeV , fixes a com¬ 
bination of the coherence length and the strength 
of the magnetic field. Approximately, these param¬ 
eters have to satisfy /c """ RhiE^), while our nu¬ 
merical calculations show that the coherence length 
should lie in the range Ic = {2 — 5) pc for acceptable 
magnetic field strengths. 

2. The shape of the knee and, for large TZmax, the 
subsequent recovery observed in the energy spectra 
of individual GR groups measured by KASGADE- 
Grande can be reproduced only, if the turbulent 
magnetic field strength is smaller than assumed e.g. 
in the JF model of Ref. [2l|, cf. Fig. |31 

3. A smaller strength of the turbulent field is also 

supported by B/G data: They constrain the gram- 
mage traversed by GRs and are consistent with a 
Kolmogorov-like power-spectrum T’{k) oc of 

the turbulent field modes. Consistency with these 
measurements at lower energies also forces us to 
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decrease i?rms- More quantitatively, we have to re¬ 
duce the turbulent field in the JF model by a factor 
~ 8, keeping the regular field unchanged. 

4. Such a reduction is in line with our determina¬ 
tion of the diffusion coefficient in a purely turbu¬ 
lent magnetic field with strength Brms = 4 /xG [i^ , 
which also disagreed by an order of magnitude with 
the extrapolation of the diffusion coefficient phe¬ 
nomenologically determined from the ratio of sec¬ 
ondary to primary nuclei. 


CR flux in the intermediate energy region up to ankle 
should be composed mainly of CRs accelerated in star- 
bust galaxies. The ankle is then interpreted as the tran¬ 
sition to another extragalactic source class, as e.g. active 
galactic nuclei or gamma-ray bursts. 

The flux of CRs escaping from starburst galaxies has 
a low-energy cutoff, when the interaction probability of 
CR nuclei on gas becomes of order one or the diffusion 
time in the intergalactic magnetic fields becomes compa¬ 
rable to the Hubble time. The magnetic horizon can be 
approximated at E'er ^ E < 10^® eV by 


Thus our model is based on relatively small values of 
the coherence length and the energy density in turbu¬ 
lent and regular magnetic fields. The first assumption 
is supported by a number of observational studies which 
derived limits on the coherence length in the Galactic 
disk of order 10 pc [s^. The second assumption ap¬ 
pears more contrived, since the required reduction factor 
/3 « 1/8 is relatively large. However, we note that the 
recent study suggests to rescale the isotropic turbu¬ 
lent field of the JF model by a factor 0.3, while it still 
predicts Faraday rotation measurements at low galactic 
latitudes that are a factor two too large. Futhermore, a 
non-uniform density of electron in the Milky Way may 
lead to an over-estimate of the turbulent Galactic mag¬ 
netic field. Additionally, one should be aware that several 
oversimplifications in our analysis may lead to a some¬ 
what too large value of /?: For instance, we have not 
properly accounted for a possible anisotropy in the tur¬ 
bulent magnetic field or a spiral structure of GR sources 
in the Milky Way. 

Let us note also that the weakness of the turbulent 
GMF in the “escape model” would have important conse¬ 
quence for the search of UHECR sources: UHECRs from 
a single source would be mainly deflected by the regular 
component of the GMF, while the spread of their arrival 
directions due to the turbulent GMF should be small. 
As a result, the search for UHEGR sources, at least in 
the case of protons or light nuclei, should be easier than 
thought before. Even for heavier nuclei, the deflections 
in the regular field of the Galaxy can be traced back in 
those patches of the sky with small turbulent fields [5l| . 
Weaker magnetic fields will also simplify the search of nu¬ 
clei sources using the methods discussed in Refs. [5^.15^. 
Thus the results of this work are an additional motiva¬ 
tion for future searches of UHECR sources, performed 
by future all-sky missions as e.g. JEM-EUSO [541 and 
KLYPVE. 0. 

B. Contribution of starbust galaxies 

It is natural to apply the ’escape model’ to other nor¬ 
mal galaxies. In particular, this model suggests that the 
CR knee in starburst galaxies is shifted by two orders 
of magnitude to higher energies [H^, since the observed 
magnetic fields of these g^axies are a factor 100 larger 
than in the Milky Way [^. Therefore, the extragalactic 
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(3) 

where Hq is the Hubble constant and the critical energy 
Ecr is defined by Rl{Eci) = Ic- 

Neither the strength nor the coherence length of the 
intergalactic magnetic field are well known. Lower limits 
on space-filling intergalactic magnetic field are around 
10“^^ [59|, while the upper limit is given by E ~ 
0.1 nG |60l| . Using for illustration B = 10“^^ G and 
Ic = Mpe, the size of the magnetic horizon at E = 10^® eV 
is for protons rhor ~ 100 Mpe and becomes comparable 
to the Hubble radius at E = 3 x 10^® eV. Since the mag¬ 
netic horizon of nuclei is smaller, a further prediction 
of this scenario is that the extragalactic flux is initially 
purely composed of protons, becoming at higher energies 
heavier and more similar to the Galactic composition. 
This is in line with our finding of an extragalactic proton 
contribution above 3 x 10^® eV, cf. Fig. [TTl and a pro¬ 
ton dominated extragalactic flux deduced from the ln(A) 
data. 

The low-energy cutoff due to interactions of CR nu¬ 
clei on gas can be estimated by rescaling Fig. [T] For 
magnetic fields a factor 100 larger, significant attenua¬ 
tion of the CR nuclei sets in below 10^® eV. Thus we 
assume that the CR flux reaching the Milky Way from 
starburst galaxies is not affected strongly by interactions 
or magnetic confinement in the for us interesting range 
E > 3 X 10^® eV. Detailed calculations of the diffuse CR 
flux from starburst galaxies will be presented in Ref. M- 


C. Connection to diffuse neutrino and 7-ray fluxes 

The recent discovery [^, of astrophysical neutri¬ 
nos with energies E > lOTeV by the IceCube experi¬ 
ment opened a new field—high energy neutrino astro¬ 
physics. The measured astrophysical neutrinos events 
are not completely uniformly distributed over the sky, 
but have an over-density towards the Galactic plane and 
a region close to the Galactic center. In Ref. [fij, it was 
suggested that these neutrinos are secondaries from CR 
interactions in the central part of our Galaxy. Later it 
was shown that the observed neutrino spectrum, which 
follows the power-law 1/E^-^® [^, has the same slope 
and normalization as the all-sky gamma-ray spectrum 
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measured by the Fermi-LAT experiment at lower en¬ 
ergies [6^. It was suggested that both spectra are dom¬ 
inated by hadronic interactions of CRs in our Galaxy. 
In this case, the CR spectrum in the central part of the 
Milky Way should be consistent with a 1/ power-law, 
which in turn agrees with the slope of the nuclei spectrum 
derived in this work. Taking into account the change of 
the power-law exponent by 1/3 in the case of Kolmogorov 
turbulence, such a spectrum is consistent with the slope 
suggested by acceleration models. 

An exception is the locally measured proton spectrum, 
which has the spectrum As argued in Ref. [^ . 

such a spectrum could by caused ’recent’ (i.e. within 
~ lOMyr) variation of the local CR flux due to a nearby 
source. Such an event might be connected to the cre¬ 
ation of the Gould belt of molecular clouds. The aged 
proton spectrum of such source would be softer than the 
average spectrum of Galactic “sea” CRs, while CR nu¬ 
clei have been spallated except of at very low rigidities. 
Finally, we note that the main contribution to the ob¬ 
served amplitude of the dipole anisotropy at E < 10^^ eV 
could be caused by this source, exceeding thereby the 
1 1E^^^ low-energy continuation of our estimate presented 
in Fig. [To] [^. 

The flux of astrophysical neutrinos contains a signif¬ 
icant fraction of neutrinos outside the Galactic plane, 
which should have an extragalactic origin. In the frame¬ 
work of the present model, these extragalactic neutrinos 
should be created by CR interaction in normal and star- 
burst Galaxies. As discussed in Ref. [^, they can ex¬ 
plain both a significant part of the diffuse gamma-ray 
background and of the IceCube signal outside the Galac¬ 
tic plane [6l|. 

D. Source rate and the slope of the injection 
spectra 

In our model, we use as average injection spectrum of 
CRs a constant power-law dN/dE ^ E~°‘' for each group 
of CR nuclei over the rigidity range 200 to 10® GV. This 
suggests that a single source class accelerates the CRs 
observed in this energy range. These sources should ei¬ 
ther all accelerate to TZmax — I x 10^^ V, or their max¬ 
imal rigidities should follow approximately a power-law 
distribution Since the required maximal energy in 
our model is high, one expects that only a subset of all 
Galactic CR accelerator is responsible for the CRs in this 
energy range. This is in line with our determination of 
the source rate, 0.4/century, that is a factor ten lower 
than the usually assumed SN rate. 

Because of the large leverage in our fits, the result¬ 
ing constraints on the exponents ai are much tighter 
than considering only, e.g., CREAM data: Rather steep 


power-laws with, e.g. ~ 2.65 for the observed pro¬ 
ton spectrum, which cannot be excluded by CREAM 
data alone are incompatible adding KASCADE and 
KASCADE-Grande data. 

VI. CONCLUSIONS 

We have shown that the knee can be entirely explained 
by energy-dependent CR leakage from the Milky Way, 
with an excellent fit to all existing data from E/Z ^ 
300 GeV to 100 PeV. In particular, all deviations from a 
single power-law behavior that are observed in the CR 
intensity of individual CR groups in the energy range 
E jZ ^ 200 GeV up to 100 PeV are consistently explained 
by rigidity-dependent CR escape. This model requires 
small coherence lengths of the turbulent field and rela¬ 
tively small turbulent magnetic fields. If these two con¬ 
ditions are fullfilled, then the CR escape time Tesc{E) ex¬ 
hibits a knee-like structure around E/Z = few x 10^® eV 
together with a recovery around E/Z 10^® eV. 

We have determined the maximal rigidity TZmax = 
Emax,Fe/(26e) ^ 10^^ V to which Galactic CR sources 
are able to accelerate CRs by identifying it with the 
energy where In (A) derived from PAO measurements 
stops to decrease. The resulting flux ratio of Galac¬ 
tic and extragalactic sources is in our model 1:1 at 
A* 2 X 10^^ eV, dropping to 0:1 at 2 x 10^® eV. The ex¬ 
tragalactic CR flux in the intermediate energy region up 
to ankle should be composed mainly of CRs accelerated 
in starbust galaxies. Since the transition from Galactic 
to extragalactic CRs happens in this model at rather low 
energies, the estimated CR dipole anisotropy is consis¬ 
tent within uncertainties with upper limits in the energy 
range 10^^ — 10^® eV, while it reproduces the measure¬ 
ments at lower energies from EAS-TOP and IceCube. 
The dipole phase is expected to change between 1 x 10^^ 
and 3 x 10^® eV, i.e. the energy range of the transition 
from Galactic to extragalactic CRs. Such a behavior cor¬ 
responds to the one observed, providing thus additional 
evidence for a transition from Galactic to extragalactic 
CRs in this energy region. 
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